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Abstract 
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A non-perturbative random-matrix theory is applied to the transmission of a 
monochromatic scalar wave through a disordered waveguide. The probability 
■ distributions of the transmittances T mn and T n = J2 m T m n °f an incident 

qq ■ mode n are calculated in the thick- waveguide limit, for broken time-reversal 

symmetry. A crossover occurs from Rayleigh or Gaussian statistics in the 
diffusive regime to lognormal statistics in the localized regime. A qualitatively 
CN ' different crossover occurs if the disordered region is replaced by a chaotic 

^ ■ cavity. 
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The statistical properties of transmission through a disordered waveguide have been 
extensively studied since 1959, when Gertsenshtein and Vasil'ev |jj computed the probability 
distribution P(T) of the transmittance T of a single-mode waveguide. It turned out to be 
remarkably difficult to extend this result to the iV-mode case. Instead of a single transmission 
amplitude t and transmittance T = |t| 2 , one then has an N x iV transmission matrix t mn 
and three types of transmittances 

N N 

T =\t I 2 T = V It I 2 T= V It I 2 (1) 

- 1 mn \ L mn\ > - 1 n / t \ L mn\ ; ± / t \ L mn\ ■ V / 

m=l n,m=l 

All three transmittances have different probability distributions, which can be measured 
in different types of experiments: If the waveguide is illuminated through a diffusor, the 
ratio of transmitted and incident power equals T/N, because the incident power is equally 
distributed among all N modes. (For electrons, T is the conductance in units of 2e 2 /h.) 
If the incident power is entirely in mode n, then the ratio of transmitted and incident 
power equals T n . For N ^> 1 this corresponds to illumination by a plane wave. Finally, 
T mn measures the speckle pattern (the fraction of the power incident in mode n which is 
transmitted into mode m). 

The complexity of the multi-mode case is due to the strong coupling of the modes by 
multiple scattering. While in the single-mode case the localization length £ is of the same 
order of magnitude as the mean free path I, the mode coupling increases £ by a factor of N. 
If iV ^> 1, a waveguide of length L can be in two distinct regimes: the diffusive regime / <C 
L ^ Nl and the localized regime L 3> Nl. The average of each of the three transmittances 
decays linearly with L in the diffusive regime and exponentially in the localized regime. In 
an important development, Nieuwenhuizen and Van Rossum |fj] (and more recently Kogan 
and Kaveh ||) succeeded in computing the probability distributions P{T mn ) and P{T n ) 
for plane-wave illumination in the diffusive regime. The former is exponential (Rayleigh's 
law) with non-exponential tails, while the latter is Gaussian with non-Gaussian tails. The 
existence of such anomalous tails has been observed in optical experiments |4j|J and in 
numerical simulations 0. From the simulations, one expects a crossover to a lognormal 
distribution on entering the localized regime. Since the theory of Refs. is based on a 
perturbation expansion in the small parameter L/Nl, it cannot describe this crossover which 
occurs when L/Nl ~ 1. 

It is the purpose of the present paper to provide a non-perturbative calculation of P(T mn ) 
and P(T n ), which is valid all the way from the diffusive into the localized regime, and which 
shows how the Rayleigh and Gaussian distributions of T mn and T n evolve into the same 
lognormal distribution as L increases beyond the localization length £ ~ Nl. We expect 
that P(T) also evolves from a Gaussian to a lognormal distribution, but our calculation 
applies only to the plane- wave transmittances T mn and T n , and not to the transmittance 
T for diffuse illumination. For technical reasons, we need to assume that time-reversal 
symmetry is broken by some magneto-optical effect. Similar results are expected in the 
presence of time-reversal symmetry, but then a non-perturbative calculation becomes much 
more involved. We make essential use of the quasi-one-dimensionality of the waveguide 
(length L much greater than width W) and assume weak disorder (mean free path / much 
greater than wavelength A). The localization which occurs in unbounded media when / < A 
requires a very different non-perturbative approach 0. 
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A related problem of experimental interest is the transmittance of a cavity coupled 
to two iV-mode waveguides without disorder. If the cavity has an irregular shape, it has a 
complicated "chaotic" spectrum of eigenmodes. At the end of the paper we compute P(T mn ) 
and P(T n ) for such a chaotic cavity and contrast the results with the disordered waveguide, 
which we consider first. 

Our calculation applies results from random-matrix theory for the statistics of the trans- 
mission matrix. This matrix t = uy/rv is the product of two unitary matrices u and v, 
and a matrix r = diag^, r 2 , . . .tjv) containing the transmission eigenvalues. It describes 
the transmission of electrons or electromagnetic radiation, to the extent that the effects 
of electron-electron interaction or polarization can be disregarded. The two plane-wave 
transmittances which we consider are 

^mtl = ^ U m.k U m.l V kn V l,n\/ T k T h = \ V kn\ 2r k- (2) 

k,l k 

We seek the probability distributions 

P(T mn ) = (6(T mn - N 2 T mn )}, (3a) 
P{T n ) = {8{T n -NT n ))i (3b) 

of the normalized transmittances T mn = N 2 T mn and T n = NT n . (These conventions differ by 
a factor l/L with Refs. @fj].) The brackets (• ■ ■) denote an average over the disorder. In the 
quasi-one-dimensional limit of a waveguide which is much longer than wide, the matrices u 
and v are uniformly distributed in the unitary group ||. The joint probability distribution 
of the transmission eigenvalues evolves with increasing L according to the Dorokhov-Mello- 
Pereyra-Kumar (DMPK) equation ||. The average can be performed in two steps, first over 
u and v, and then over the transmission eigenvalues T k . 

The first step was done by Kogan and Kaveh ||. The result is an expression for the 
Laplace transform of P(T n ), 

POO 

F(s)= dT n exp(-sT n )P(T n ), (4) 
Jo 

which in the thick-waveguide limit (N — > oo, L/l — > oo, fixed Nl/L) is exactly given by 

F(s) = (u k (l + sr k )- 1 ). (5) 

The same function F(s) also determines P(T mn ), which in the same limit is related to P{T n ) 
by i 

P(T mn ) = | o °° dT n T- 1 exp (-T wn /T n ) P(TJ. (6) 

The next step, which is the most difficult one, is to average over the transmission eigen- 
values in Eq. (^j. The result depends on whether time-reversal symmetry is present or not 
(indicated by (3 = 1 or 2, respectively). In Refs. @|§, InF was evaluated to leading order in 
L/Nl, under the assumption that the waveguide length L is much less than the localization 
length £ ~ Nl. Here we relax this assumption. 
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exactly. This was shown by Rejaei [IT 



super symmetric non-linear a model [12 



We consider the case of broken time-reversal symmetry {(3 = 2). Then the probability 
distribution of the tj-'s is known exactly, in the form of a determinant of Legendre functions 
P v (TI| . Still, to compute expectation values with this distribution is in general a formidable 
problem. It is a lucky coincidence that the average (|5|) which we need can be evaluated 

using a field-theoretic approach which leads to a 
It was recently proven Jl3| that this super sym- 
metric theory is equivalent to the DMPK-equation used in Ref. From Rejaei's general 
expressions we find 

F(s) = l-2sJ2 / dk f p (k) tanh(lnk)P Hik _ 1) {l + 2s), 
(2p + l)k 



f P (k) 



cxp 



L[(2p+l) 2 + k 2 



(2p + l) 2 + k 2 ~' r \ ANl 
Inversion of the Laplace transform (|J) yields P(T n ), 

o 2K hik ar n 



^ poo 

P{Tn) = Y. / dkf p {k)sm\y{\-Kk) 



dT n (7r s T n e^y/r 



(7) 



where K v is the Macdonald function. One further integration gives P(T mn ), in view of Eq. 
(§). Results are plotted in Fig. [I]. The large T n and T mn tails are 



P{T m 
P{%, 



2 V Tm 



T mn > 1, (Nl/Lf 
% > l,Nl/L. 



(8a) 
(8b) 



It is worth noting that Fyodorov and Mirlin |T|| found the same tail as Eq. (|8aD for the 
distribution of the local density of electronic states in a closed disordered wire. It is not 
clear to us whether this coincidence is accidental. 

The diffusive and localized limits can be computed from Eq. (|5|) by using the known 
asymptotic form of the distribution of the Tfe's. In contrast to the full result ([?]), which holds 
for j3 = 2 only, the following asymptotic expressions hold for any j3. In the diffusive regime, 
for L Nl, we may expand InF in cumulants of the linear statistic A = X^ m (l + st^): 



InF(s) = ln(e 



(A) + |Var A + 0(L/NI). 



(9) 



The mean and variance of A can be computed from the general formulas of Refs. [|IT)|,|il^,|To| : 



(A) — — asmlrVs + 



Var A 



1 

1 




(10a) 
(10b) 



valid up to corrections of order L/Nl. To leading order in L/Nl one has the /3-independent re- 
sult of Refs. yielding Gaussian and Rayleigh statistics for L/Nl — > 0. The /3-dependent 
terms in Eqs. ([II]) are the first corrections due to localization effects. In Fig. ^|we plot P{T n ) 
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resulting from Eqs. @ and fllCf) . The /3-independent result of Refs. pm (not shown) is very- 
close to the (3 = 2 curve. This figure indicates that the /3-dependence is essentially quanti- 
tative rather than qualitative. 

In the opposite, localized regime (L ^> Nl), only a single transmission eigenvalue con- 
tributes significantly to Eq. (H). This largest eigenvalue r has the lognormal distribution 



\8ttL J 



8L \pm 



(11) 



It follows that lnT mn and In T n are also distributed according to Eq. (JTTj) in the localized 
regime. The approach to a common lognormal distribution as L/Nl increases is illustrated 
in Fig. |, using the exact (3 = 2 result of Eq. (0). 

We contrast these results for a disordered waveguide with those for a chaotic cavity, 
attached to two iV-mode leads without disorder. Following Ref. |TJ| we assume that the 



2N x 2N scattering matrix of the cavity is distributed uniformly in the unitary group if 
(3 = 2 or in the subset of unitary and symmetric matrices if /3 = 1. Then P(T mn ) and P{T n ) 
follow from general formulas [0 for the distribution of matrix elements in these socalled 
"circular" ensembles. For (3 = 2 the result is 



P(T mn ) = (2N - 1) (1 - T mn ) 2N ~ 2 , (12a) 
P(T n ) = W 2 5) PT»(1 " Tn)]"- 1 • (12b) 



For (3 = 1 Eq. QT2aJ) should be multiplied by ±F(N - |,1;2JV-1;1- T mn ) and Eq. 



by \F[N — ~, 1; iV; 1 — T n ), where F is the hypergeometric function. These are exact results 
for any iV. If iV — > oo, P(T mn ) is an exponential distribution with mean 1/2N, and P(T n ) is 
a Gaussian with mean 1/2 and variance 1/8N. This is similar to the disordered waveguide, 
with iV playing the role of Nl/L. As shown in Fig. |], the distributions for N of order unity 
are quite different from those in a disordered waveguide with Nl/L of order unity. For N = 1 



the distinction between T mn ,T n , and T disappears and we recover the results of Ref. fig . 

In conclusion, we have presented a non-perturbative calculation of the distributions of 
the plane-wave transmittances T mn and T n through a disordered waveguide without time- 
reversal symmetry, which shows how the distributions cross over from Rayleigh and Gaussian 
statistics in the diffusive regime, to a common lognormal distribution in the localized regime. 
Qualitatively different distributions are obtained if the disordered region is replaced by a 
chaotic cavity. Existing experiments have been mainly in the regime L Nl where the 
perturbation theory of Refs. [|2]|| applies. If the absorption of light in the waveguide can be 
reduced sufficiently, it should be possible to enter the regime L ~ Nl where perturbation 
theory breaks down and the crossover to lognormal statistics is expected. 

This research was supported by the "Nederlandse organisatie voor Wetenschappelijk On- 
derzoek" (NWO) and by the "Stichting voor Fundamenteel Onderzoek der Materie" (FOM). 



5 



REFERENCES 



[1] M. E. Gertsenshtein and V. B. Vasil'ev, Teor. Veroyatn. Primen. 4, 424 (1959); 5, 3(E) 

(1960) [Theor. Probab. Appl. 4, 391 (1959); 5, 340(E) (I960)]. 
[2] Th. M. Nieuwenhuizen and M. C. W. van Rossum, Phys. Rev. Lett. 74, 2674 (1995). 
[3] E. Kogan and M. Kaveh, Phys. Rev. B 52, 3813 (1995). 
[4] A. Z. Genack and N. Garcia, Europhys. Lett. 21, 753 (1993). 

[5] J. F. de Boer, M. C. W. van Rossum, M. P. van Albada, Th. M. Nieuwenhuizen, and 

A. Lagendijk, Phys. Rev. Lett. 73, 2567 (1994). 
[6] I. Edrei, M. Kaveh, and B. Shapiro, Phys. Rev. Lett. 62, 2120 (1989). 
[7] B. L. Al'tshuler, V. E. Kravtsov, and I. V. Lerner, in Mesoscopic Phenomena in Solids, 

edited by B. L. Al'tshuler, P. A. Lee, and R. A. Webb (North-Holland, Amsterdam, 

1991); B. A. Muzykantskii and D. E. Khmel'nitskh, Phys. Rev. B 51, 5480 (1995); 

V. I. Fal'ko and K. B. Efetov, preprints (|cond-mat/ 95030961 , 9507091); A. D. Mirlin, 

preprints frond-mat/ 95U7U84 9508093). 
[8] P. A. Mello, E. Akkermans, and B. Shapiro, Phys. Rev. Lett. 61, 459 (1988). 
[9] O. N. Dorokhov, Pis'ma Zh. Eksp. Teor. Fiz. 36, 259 (1982) [JETP Lett. 36, 318 (1982)]; 

P. A. Mello, P. Pereyra, and N. Kumar, Ann. Phys. (N.Y.) 181, 290 (1988). 
[10] C. W. J. Beenakker and B. Rejaei, Phys. Rev. Lett. 71, 3689 (1993); Phys. Rev. B 49, 

7499 (1994). 
[11] B. Rejaei, Phys. Rev. B (to be published). 
[12] K. B. Efetov, Adv. Phys. 32, 53 (1983). 

[13] P. W. Brouwer and K. Frahm, preprint (|cond-mat/95080"56| ) . 

[14] Y. V. Fyodorov and A. D. Mirlin, Int. J. Mod. Phys. B 27, 3795 (1994). 

[15] C. W. J. Beenakker, Phys. Rev. B 49, 2205 (1994). 

[16] J. T. Chalker and A. M. S. Macedo, Phys. Rev. Lett. 71, 3693 (1993); Phys. Rev. B 
49, 4695 (1994). 

[17] A. D. Stone, P. A. Mello, K. A. Muttalib, and J.-L. Pichard, in Mesoscopic Phenom- 
ena in Solids, edited by B. L. Al'tshuler, P. A. Lee, and R. A. Webb (North-Holland, 
Amsterdam, 1991). 

[18] H. U. Baranger and P. A. Mello, Phys. Rev. Lett. 73, 142 (1994); R. A. Jalabert, J.-L. 

Pichard, and C. W. J. Beenakker, Europhys. Lett. 27, 255 (1994). 
[19] P. Pereyra and P. A. Mello, J. Phys. A 16, 237 (1983). 



6 



FIGURES 



FIG. 1. Distributions of (a) T ri = NT n for L/Nl = 0.05, 0.1, 0.5, 1.5, 2.0, and 2.5, and (b) 
T mn = N 2 T mn for L/Nl = 0.05, 0.5, 2.5, 5, and 10. Computed from the exact j3 = 2 expressions 
(|6|) and (0). The dotted curves are the limits L/Nl — > of an infinitely narrow Gaussian in (a) 
and an exponential distribution in (b) (note the logarithmic scale). The inset in (b) shows the 
waveguide geometry considered (disordered region is shaded). 



FIG. 2. Distribution of T n calculated from the perturbation expansion (Q), ()10|), for (3 = 1,2 
and L/Nl = 0.1,0.5. 

FIG. 3. Distributions of T n and T mn for (3 = 2 and L/Nl = 5, 10, 20, computed from Eqs. @ 
and (0). The dotted curve is the lognormal distribution ([ll]) which is approached as L/Nl — > oo. 

FIG. 4. Distribution of T„, for a chaotic cavity attached to two A^-mode leads (inset). The 



curves are computed from Eq. (12b), for j3 = 1, 2 and N = 1, 2, 20 
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Figure la. 




Figure lb . 
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Figure 4. 



